A Condensation of Interacting Bosons in Two Dimensional Space 



S. G. Rajeev 
Department of Physics and Astronomy 
University of Rochester 
Rochester,NY 14627 

Abstract 

We develop a theory of non-relativistic bosons in two spatial dimensions with 
a weak short range attractive interaction. In the limit as the range of the 
interaction becomes small, there is an ultra-violet divergence in the problem. 
We devise a scheme to remove this divergence and produce a completely finite 
formulation of the theory. This involves reformulating the dynamics in terms 
of a new operator whose eigenvalues give the logarithm of the energy levels. 
Then, a mean field theory is developed which allows us to describe the limit 
of a large number of bosons. The ground state is a new kind of condensate 
(soliton) of bosons that breaks translation invariance spontaneously. The 
ground state energy is negative and its magnitude grows exponentially with 
the number of particles, rather than like a power law as for conventional many 
body systems. 
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There is a well-known theory of Bose gases in one dimensional space with 
a delta- function interaction However the two-dimensional bose system 

with attractive delta-function interactions is logarithmically divergent even 
in the case of two bodies ||,^],||,f?]]- The hamiltonian of the theory is scale 
invariant so that the ground state energy is either zero or negative infinity; 
in the case of attractive interactions it is negative infinity. The two-body 
problem has been intensely studied even beyond perturbation theory and it 
is now understood how this divergence can be removed by a renormalization. 
Recently we extended this analysis to the case of the three body problem; 
while not exactly solvable, the problem can be reformulated in a finite form 
without using perturbation theory 0. 

In this paper (for a more detailed discussion see Ref. ||) we will describe 
a non-perturbative, fully renormalized formalism for the corresponding non- 
relativistic field theory; i.e., the case of an arbitrary number of bosons. We will 
then develop a mean field theory to obtain the ground state of the system. 
The results are striking: the ground state is a new kind of condensate (or 
soliton, a bound state of a large number of particles) of bosons, which breaks 
translation invariance spontaneously. Moreover, the binding energy grows 
exponentially with the number of particles. 

Of course, such a study of the binding of bosons is only possible because 
we now have a non-perturbative renormalization method for this theory. This 
bosonic field theory is asymptotically free (i.e., the running coupling constant 
decreases logarithmically with energy) and hence provides a good testing for 
our new ideas of renormalization |@j. These ideas apply also to other asymp- 
totically free theories such as the two-dimensional non-abelian Thirring model 
and perhaps, even to four dimensional gauge theories such as Quantum 
ChromoDynamics . 

We begin with the naive hamiltonian of the system of bosons with an 
attractive self- interaction of zero range p| : 



H = \ [ \V0\ 2 d 2 x-g [ ft 2 (x)j> 2 (x)d 2 x = [ £^t(p)0(p) 



dp 




dpidp2dp' 1 dp' 2 



(27T) 2 6( Pl +p 2 - V \ -P' 2 )^{PIW{P2)4>{P I 1 )^P'2)- 



(2tt) 8 



(1) 



Here, (j)(x) is a complex scalar field in two dimensional space satisfying the 



usual canonical commutation relations; (j)(p) = f e tp ' x (p(x)dx is its Fourier 



transform. In the units we are using, (j) has dimensions of (length)" 1 and 
H has dimensions of (length) -2 ; it follows that the coupling constant g is 
dimensionless. This means that the lowest eigenvalue of H (in a sector with 
fixed number of particles ) is either zero or infinite. If g > 0, corresponding to 
an attractive potential, it is known that the ground state energy is — oo even 
in the case of two particles §. 

Following the standard ideas of field theory, we must regularize the prob- 
lem by, for example, introducing a cut-off in the momentum; also, the coupling 
constant must be made a function of the cut-off in such a way that the ground 
state energy is finite. This is not always possible; for example the analogous 
problem in three dimensions requires an infinite number of coupling constants 
to be adjusted this way. When it is possible to remove all the divergences, by 
making a finite number of coupling constants depend on the cut-off, we say 
that the theory is renormalizable. We will show that our bosonic field theory 
in two dimensions is renormalizable. Moreover, we will produce a completely 
finite reformulation of the theory without the use of perturbation methods. 
This allows us to ask questions about the bound states of the system, even 
the bound state of a large number of particles, a more detailed account of the 
methods used is in Ref. ||. 

The regularized hamiltonian is, in momentum space, 



H A = H + H lA 



where 




and 



u ( \ \ f dpidp 2 dp' 1 dp' 2 , , 

#1A = -5(A) J T^jTg PA{Pl ~ P2)PA(Pi ~ P 2 ) 

(2vr) 2 5( Pl +p 2 - p[ - p' 2 )^(pi)4>Hp2)<P(Pi)4>(P2)- 

Here p\ is a cut-off function that is equal to one near the origin in momentum 
space and vanishes rapidly at infinity. For example we could choose 

pa(p) = e(\ P \ < A). 

We will see that this cut-off in the relative momentum is sufficient to regularize 
the theory. 

The solution of a quantum mechanical system is equivalent to determining 
the resolvent (Green's function) of its hamiltonian. Thus we would have an 
exact renormalization of our theory if we can construct the resolvent (H& — 
E)^ 1 in the limit as A — > oo. Our key result is that such a formula for the 
resolvent can be obtained, 

lim — -1— := R(E) 

A— >oo H\ — hi 

i i i lt r i , -e 1-1 
+ -— — -w 



log — + w 



2tt p 



Hn — E 



Hq — E 2 Hq — E 

We will give below explicit formulas for the operators b and W. This gives 
an indirect construction of the hamiltonian of the renormalized theory as the 
operator for which R(E) = [H — E] -1 . We cannot get an explicit formula for 
the hamiltonian itself; instead we express its resolvent in terms of the resolvent 
of an explicitly known operator, W. The dynamical information about the 
system is thus encoded in the integral operator W ('the principal operator'); 
it plays as important a role in our theory as the hamiltonian does in usual 
physical theories. Once we reformulate the theory this way, the standard 
approximation methods of quantum theory can be applied to the principal 
operator. Later on we will in fact apply mean field methods to get the ground 
state energy of the many body system. 

The operators Hq, b, W involve no dimensional parameters; the only 
parameter in the theory is p. For example the bound state energy levels of 
the system are given by the eigenvalues w a of W: 



E a = -[i 2 e~ 2 ™«. (7) 

Once we determine one eigenvalue, it fixes \i and then all other eigenvalues 
are determined. For example, — fj? is the ground state energy of the two- 
body system; once we know that the ground state of the n-body system for 
any n is determined. The original dimensionless coupling constant has been 
traded for the dimensional constant \x in the process of renormalization: this 
is dimensional transmutation |l0| ]. 

In order to get the above formula for the resolvent, we have to overcome 
some technical difficulties. The coupling constant appears in H\ multiplica- 
tively, which makes it difficult to perform a renormalization except, of course, 
in perturbation theory. We will introduce some auxilliary variables that will 
help us to rewrite the the theory in such a way that coupling constant ap- 
pears additively rather multiplicatively. This is key to our method for exact 
renormalization. Define creation-annihilation operators x{p)jX.Hp) satisfying 
the product (not commutation relations) 

x(p)xV) = (27r) 2 5(p-p'), x(p)x{p') = = X ] (p)x ] (p')- (8) 

In addition, these operators commute with the bosonic fields. These rela- 
tions are the q — > limit of the g-oscillators studied by many authors (e.g., 
PH ). These relations have an obvious representation on C © L 2 (R 2 ): there 
is an empty state |0 > and states containing one auxilliary particle such as 
J v(p)xHp) (2^)2 10 > - Since the product of the creation operators vanish, 
there can no more than one such particle in any state: far more restrictive 
than the Pauli exclusion principle for fermions. We dont have to attach a 
physical significance to these new auxilliary particles: they can be viewed as 
a mathematical device that simplifies the analysis. [] 



1 However there are examples (quarks or ghosts) of such auxilliary variables that are important 
physically as well. We propose that the auxilliary particles created by our operators \i x be called 
angels. 
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On the combined Hilbert space B' of bosons and the auxiliary particles 
(B' = B © B <8> L 2 (R 2 ), where B) is the usual bosonic Fock space), we define 
the operator: 



H' A = H U + J [dp 1 dp 2 dp 3 ]pA(pi - Pi) 

4> ] {Pi)4> [ {P2)x{Pz){^) d KPi +P2 - Pa) + h.c. 



+ 



5(A) 



ni. 



Here, IIo, III are the projection operators to the subspaces containing zero or 
one auxilliary particle respectively: 



n = J [dp]x(p)xHp), ni = J [dp]xHp)x(p)- 

We are interested in getting a formula for the resolvent of the hamiltonian 
H\. We will find it convenient to think in terms of the auxilliary hamiltonian 
H' A and and the related operator R' A = (H' A — EHq)^ 1 . 

Let us split the operators into 2x2 blocks according to auxilliary particle 
number: 

fa &t\ 1 fa ft\ 

H' A - EU = , R' A {E) = — - = 

V b d J H A -EU \p 5 J 



with 



a, a 



:B^B, b\ft :B®L 2 {R d ) ->■ B,d,5 : B ® L 2 (R d ) B <g> L 2 (R d ). 



Two different ways of writing the elements of the inverse matrix will be useful 
(the proof is elementary) : 



_ bU^b] 1 = a = a' 1 + crV[d - ba^b^ba,- 1 , 



-i 



d - ba- 1 ^ =5 = d- 1 + d-H[a - bU^b^bU' 1 , 



- 5ba~ l = /3 = -d~ l ba. 
Using the left hand side of the formula for a we can see that 

a = (H A - E)- 1 
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is the just the resolvent of the original bosonic system. The other way to 
write this quantity gives us 



+ 



1 



H\ — E Hq — E 2 Hq — E 



-tf$ A (E)- l b 



Hq — E 



Here 
and 



1 



5(A) 
xHpi +P2 



dpidp2dp' 1 dp' 2 
(2^)8 

Hpi)<P(p2) 



Ph{Pl -P2)pk{Pl -P' 2 ) 
1 



Hq — E 



^{P\)^{P2. 



X(Pl +P2)- 



6 f = J [dpidp 2 }p\(pi - P2)<t>\pi)4>\p2)x{Pi +P2)- 

The point of this formula is that it expresses the resolvent of the hamil- 
tonian in terms of that of the free theory and the operator &\(E). Moreover, 
the coupling constant appears additively in <3>a; we will be able to separate 
the divergence in <I>a and remove it by choosing g(A) to be an appropriate 
function of A. Then the limiting operator $>(E) = liniA^oo ®a{E) will exist. 
The roots of the equation &(E)\s >= will give the energy levels. Thus we 
will obtain a completely finite reformulation of the problem: what we have 
called elsewhere a transfinite formulation. It is important that no where in 
this argument we will use perturbation theory. 

In more detail, we can use the canonical commutation relations to normal 
order the operators in <E>a- This will give 

§a(E) = cT^A) - J [dp 1 d P2 dp 1 dp' 2 }p A (p 1 - p 2 )pa(p'i - P 2 ) 

X\P1 +P2) 

1 

-<p{pi)(t){p 2 ) 



+4(2n) d S(p 1 -p[)ct>\p 2 )- ' 



Hq + w(pi) + U(jp' 2 ) + L0( P2 ) - E 



<P(P2 



+2{2i:) d 5(p 1 -p' l ){2^ d 5(p 2 -p' 2/ 

x(p'i+p' 2 )- 



flo + w(pi) +oj{p' 2 ) -E 



The only divergent term is the first one the square brackets. The divergence 
is removed if we choose 



p 2 + /j, z 



Here we are trading the dimensionless coupling constant g for the dimensional 
constant /i. The renormalized theory will depend on fi alone and not the cut- 
off A or on the coupling constant: dimensional transmutation. 
The operator &\(E) has a limit as A — > oo: 

,2 J2 



$(E) = J [dp] X Hp)^,H + P -- E)x(p) 
-J [dp 1 dp 2 dp' 1 dp' 2 ]x l (pi +P2) 

i 

$ {pi)4>Hv2)- 



H + lo(p[) + u(p' 2 ) + ufa) + lu( P2 ) - E 
+A{2-K) d 5( Px -p' 1 )<t>\p' 2 ) 1 



XiPl +P2)- 



H + u(p[) + u(p' 2 ) + lo{ P2 ) - E 



<j)(pi)<j>(p 2 



(21) 



(22) 



Here, 



?l 2 ' 2 ; 



1 ^ 
-; — m — ^ 

47T fl 2 



P 2 + fi 2 p 2 + Z/ 2 



is a convergent integral. We can rescale all the momemta by y/\E\ to get 
$(p,E) 



-!-ln-^ + ^ 



(From this point on our momentum variables p, p' and position variables x 

2 

are dimensionless.) The principal operator W is, (cj(p) = \) 



ir = 27r ' ^ xt ^ log 



^o+u;(p) + l 
[dpidp'zdp^dp'^xHpi +P2 



xip) 



H + lo(p[) + w(p£) + w(pi) + w(P2) + 1 



+4(2 7 r)^(p 1 -y i )^(P 2 ) 
x(pi + P 2 )- 



ff + w(pi) + w(p 2 ) + w(p 2 ) + 1 



0(jPl)0(P2 
0(P2) 



(23) 



(24) 



(25) 



2 We assume that £7 < so that it is not in the spectrum of Hq. The resolvent is defined elsewhere 
by analytic continuation. 
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This proves the formula @ for the resolvent. 

Thus the eigenvalues of the operator W will determine the energy levels. 
In fact, 

W\tp >= w\ip >, E = -n 2 e~ 2ww . 

The problem is completely well-posed in terms of the operator W; it is a re- 
placement for the hamiltonian of the theory. In fact there are many subtleties 
in the definition of the hamiltonian itself (e.g., domain self-adjointness), all 
of which can be avoided by working with W. We call this remarkable new 
operator that embodies the dynamics of the system the 'Principal Operator'. 
Roughly speaking it represents the logarithm of energy. 

Now that we have a finite formulation of the theory, we can look for 
analogues of all the approximation methods of quantum theory. In particular, 
we can look for a mean field theory. The idea is to use a variational principle 
to estimate the lowest eigenvalue of W in the sector with n bosons. Now 
W acts in the Hilbert space of n bosons and one auxilliary particle. The 
variational ansatz corresponding to mean field theory is the one where all the 
bosons are in the state u and the angel is in the state v: 
dpi--- dp n dp 



It, V >- 



-u(pi) ■ ■ ■ u(p n )v(p)4> ] (pi) • • • t (p n )X t (p)|O > • 



( 2vr )2n+2 

The expectation value of W in this state becomes (for large n) the 'principal 
function' 



U 



1 

2^ 



where 



v(p)\ 2 log[nho(u) + u)(p) + 1] 

[dpidp 2 dp' 1 dp l 2 }v{pi +P2)*v(p' 1 +p 2 ) 
u*{p' 1 )u*{p' 2 )u(pi)u(p2)n(n - 1) 
nh (u) + uj(p[) + u(p' 2 ) + u(pi) + Ld(p 2 ) + 1 
(2tt) 2 5(pi - p'i)nu*{p' 2 )u(p 2 ) 
nh (u) + uj(p[) + U)(p' 2 ) + oo(p 2 ) + 1 



ho(u) = / \u(p)\ 2 u;(p)[dp]. 



+4 



This U is to be minimized subject to the normalization conditions 



/ 



u{p)\ 2 [dp] 



\v(p)\ 2 [dp] = 1. 



If we keep just the leading terms as n — > oo, 

2 



T 



-n 



Jv*(pi + P2)u{pi)u{p 2 )[dp 1 dp 2 



h (u) 



+ — logn-— log/i(«) + 0(n x ). 



The functional (which is independent of n) can also determined as above, 
but it takes a bit more work. Thus we have the approximate variatonal 
principle for the ground state energy of a system of n bosons: 



where 



27TT1 

E n = -f—[C 1 + 0{-)] 
n n 



J>(p) 1 2 [dp] IHp)\ 2 [dp] J u{p)\u{p)\ 2 [dp] 



f = inf 

u,v 

Jv*(pt + P2)u(pi)u(p 2 )[dp 1 dp 2 } 

Also, C\ is a constant; it is not necessary to determine it for the leading 
behavior in the large n limit of energy. We can eliminate v rather easily to 
get an equivalent form (written in terms of the position space wavefunction 
fi(*)=/u(p)e*»^)): 

£ = MI[u] 



where, 



J \u(x)\ 2 dx J l\Vu(x)\ 2 dx 



j\u{x)\ A dx 

All the variables in our problem are dimensionless. The minimization of 
log /[it] gives the partial differential equation 



V 2 u- pu + X\u\' z u = 



where, 



P 



J\Vu\ 2 d 2 x 



X 



J\Vu\ 2 d 2 x 



J \u\ 2 d 2 x ' J \u\' i d 2 x 

This nonlinear Schrodinger equation has been shown to have solutions in 
[12], Theorem 6.7.25. ( Because of the scale covariance it is enough to find 



10 



solutions with one positive value of (3 and A.) We expect the ground state to 
have at least circular symmetry; we get a sort of nonlinear Bessel's equation: 

with u{x) = y/rv(r) and r = |x|. It is easy to solve this equation numerically; 
we plot the result in the figure. The value of £ we get is about 12. The 
solution of course breaks translation invariance spontaneously but not rotation 
invariance. It is in fact a 'soliton' in the sense of field theory. 

Ground State Wave Function 

4 




2 4 6 8 

It follows that the ground state energy (binding energy) grows exponen- 
tially with the number of particles: 

Here, — /U 2 is the binding energy of the two boson system; the dots denote 
terms of lower order in the number of bosons. At zero temperature, essentially 
all the available particles will fall into this bound state. That is the sense in 
which this state is a condensate. 

The wavefunction was determined above in terms of a dimensionless po- 
sition variable. By putting back the factor y/\E\ we scaled out, we see that 
the size of the condensate behaves like \E n \~~z ~ [i~ x e~^ n \ it shrinks with 
increasing number of particles. Although the results above are asymptotic for 
large n, we expect (by analogy with mean field theory in atomic physics) that 
they will hold true even for moderate values of n of about 3 or more. 
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It may be possible to observe such a condensation in cold optical traps 
of dilute bose gases, if we can trap the atoms close to some surface p3| . 
Even without a confining potential in the plane of the surface, the atoms will 
form a 'self-trapped' bound state. It may be possible to tune the attractive 
force between atoms using Feshbach resonances |14| . Then, even with an 
exponentially small size for the condensate, the atoms can be far apart and 
the approximation of zero range will be valid. Whether large particle numbers 
in the condensate can be achieved is not clear yet, however. There might also 
be a realization in terms of vortices in superconductors which are on the 
boundary between type I and type II |8). 

In the case of the three dimensional attractive Bose gas, even after the 
above renormalization the energy is not bounded below ||; our methods fail. 
In fact there are indications ( both theoretical and experimental) that in this 
case, the system is truly unstable [15,19]. 
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